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Abstract. We use a unified analytical formulation - de- 
veloped in previous papers - for the multiplicity functions 
of clusters and galaxies. This method is free from the 
cloud-in-cloud problem encountered in earlier approaches 
and well adapted to the description of the non-linear clus- 
tering features. It is especially suited to simultaneously de- 
scribe rich clusters, groups and galaxies, consistently with 
the hierarchical picture of gravitational clustering, as well 
as their evolution in time. We find a good agreement with 
observations for the cluster temperature function and we 
compare our method with the standard Press-Schechter 
prescription. We also obtain the main properties of the 
Sunyaev-Zel'dovich effect (mean and variance). Then, us- 
ing a simple model for the cluster X-ray luminosity (taking 
into account entropy considerations), we obtain the X-ray 
luminosity distribution of groups and clusters. 

Then, using the same formalism we derive the galaxy 
and quasar multiplicity functions. In particular, we show 
that the use of the standard Press-Schechter prescription 
leads to erroneous conclusions at low redshifts while our 
approach provides a reasonable agreement with observa- 
tions in a natural fashion because it is able to distin- 
guish galactic halos from groups or clusters. Finally, we 
derive the contribution of quasars and galaxies to the X- 
ray counts. Thus, we obtain a global and consistent picture 
of the X-ray emissions from all structures. In particular, 
we show that future observations (e.g., from AXAF) could 
provide interesting information on galaxy evolution. In- 
deed, they will constrain the importance of a possible hot 
diffuse gaseous phase in galactic halos and they could re- 
veal massive galaxies which are just being formed, through 
the X-ray emission of their cooling gas. 



Key vifords: galaxies: clusters - galaxies: evolution - 
galaxies: mass function - quasars: general - cosmology: 
large-scale structure of Universe - cosmic microwave back- 
ground 



1. Introduction 

Clusters of galaxies are the largest virialized objects, char- 
acterized by scales which have just entered the non-linear 
regime. Thus, they are very rare and their number den- 
sity is very sensitive to the cosmological parameters and 
to the amplitude of the initial density fluctuations (within 
the framework of the usual hierarchical scenarios). As a 
consequence, as proposed by Oukbir & Blanchard (1992) 
many authors have compared observations with predic- 
tions from numerical simulations (e.g.. Eke et al. 1996) or 
analytic approaches (e.g., Oukbir & Blanchard 1997) in 
order to obtain constraints on (ilm, ^^a) and on the initial 
power-spectrum P{k). However, the cluster temperature 
- X-ray luminosity relation involves non-gravitational ef- 
fects since simple scaling laws (Kaiser 1986) recovered by 
numerical simulations (Eke et al. 1998) that neglect ra- 
diative cooling and supernova or quasar feedback disagree 
with observations. This could be due to a preheating of the 
IGM by QSOs or supernovae (Valageas & Silk 1999b; Cav- 
aliere et al. 1997; Ponman et al. 1999) which has not been 
considered yet by numerical simulations. Thus, analytic 
approaches are still needed in order to describe clusters. 
Moreover, they explicitly show the connection of cluster 
characteristics (e.g., their mass function) with other fea- 
tures of the universe (e.g., various properties of the under- 
lying density field, of galaxies or quasars). 

The standard way to estimate the mass multiplicity 
function of collapsed objects by analytical means, is the 
well-known Press-Schechter (1974) approximation - here- 
after PS - that directly recognizes in the initial, linear, 
fluctuations of the density field the overdensities that will 
eventually collapse. This approximation, however, suffers 
from many drawbacks. One is that it is customary -and 
necessary in order to bring the analytical form in agree- 
ment with numerical simulations- to multiply the result by 
an ad-hoc factor of two. This cannot be justified by the 
standard excursion sets argument for realistic filters like 
the top-hat in real-space as shown in Valageas & Schaef- 
fer (1997) - hereafter VS - (see also Peacock & Heavens 
1990). Indeed, in such a case the excursion sets imply a 
renormalization factor which goes to unity at large masses. 
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Another fundamental problem, inherently related to 
hierarchical clustering, is to describe objects embedded 
within other, less dense but nevertheless virialized, ob- 
jects, that is to describe subclustering. This has long been 
known (Bardeen et al. 1986) not to be reliable within ap- 
proaches directly based on recognition of the linear over- 
densities that will form objects, and is called the cloud- 
in-cloud problem. Indeed, an approach based on count- 
ing overdensities in the linear regime leads to severe over- 
counting, so severe that the same objects are erroneously 
counted an infinite number of times (VS). Thus, to de- 
scribe a universe made of dense galaxies (10"* times the 
mean density), as well as of galaxy clusters (10^ times 
the mean density), is out of reach of such approaches, 
while the description of intermediate objects such as Lya 
clouds, that may even have densities below the mean, is 
unthinkable within this framework. Indeed, note that the 
PS approximation can only deal with just-virialized halos. 

Recent progress (VS) in the description of the non- 
linear density field, through a non-linear scaling model 
based on earlier approaches (Schaeffer 1984, 1985; Balian 
& Schaeffer 1989; Bernardeau & Schaeffer 1991), with an 
understanding of its relation to the initial spectrum of fluc- 
tuations, leads to analytical expressions that enable one 
to directly count the overdensities in the actual non-linear 
density field. The result exhibits many analogies with the 
PS approximation, but it gives a correct normalization 
and it also solves the cloud-in-cloud problem (there are 
no divergences) . It also provides the correlations of these 
objects and the associated bias as compared to the matter 
distribution (Valageas et al. 2000b; Bernardeau & Schaef- 
fer 1992,1999). 

For the purposes of the present paper, the scaling 
approach gives a clear answer to the following problem. 
Choose any density contrast A = {p — p)/p and define 
objects as having an inner density contrast larger than 
this threshold. Then all the mass in the universe (except 
a nonvanishing but negligible fraction in the very under- 
dense regions) lies within such objects, that have a dis- 
tribution of mass or size given by the theory. The choice 
of a larger contrast results in a different partition of the 
same mass into smaller objects, the procedure being valid 
for arbitrarily large contrasts (with obvious limits at the 
kpc scale where pressure and angular momentum come 
into play) as long as one remains in the non-linear regime. 
This is the way subclustering can be properly accounted 
for. Then, the same theory can simultaneously describe 
virialized clusters, which have a contrast of ~ 200, as well 
as galaxies (Valageas & Schaeffer 1999), whose contrast is 
~ 5000, embedded or not in the latter objects. The proce- 
dure also holds for small density contrasts, even negative, 
provided one is at scales (in the present universe, for in- 
stance, not much above 1 Mpc) where the correlation func- 
tion is large enough so as to insure one is in the fully non- 
linear regime. This corresponds to underdense non-linear 
objects embedded in much less dense, nearly void regions. 



The new feature introduced by the scaling model is that, 
because one directly works in the non-linear regime, this 
separation can be done properly. The same procedure ap- 
plied to future non-linear objects in the linear regime as 
is done in the PS approach would lead to the above men- 
tioned divergences (VS), which is another indication that 
defining objects directly in the linear regime, as is done in 
the latter approach, is unsecure. 

The analytical predictions of the scaling model have 
been checked against numerical simulations (Valageas, 
Lacey & Schaeffer 2000a) under the most extreme con- 
ditions, searching for non-linear objects with density con- 
trasts ranging from values as large as A = 5000 down to 
negative values A = —0.5, spanning four orders of mag- 
nitude and limited only by the accuracy of the simula- 
tions. Thus, this approach provides reasonable results, at 
least in the above range, allowing us to describe non linear 
objects and their evolution. These theoretical predictions 
are also in the line of the findings of Moore et al. (1999a) 
that resolve some of the structure of dark matter halos 
in their simulation. As a consequence, the mass functions 
provided by this approach give, with the same parame- 
ters, a unified description of very different objects, such 
as galaxies and quasars (Valageas & Schaeffer 1999), Lya 
clouds (Valageas, Schaeffer & Silk 1999), the ionization 
flux emitted by QSOs (Valageas & Silk 1999a), allowing 
for a consistent picture of the reheating and the reioniza- 
tion history of the universe (Valageas & Silk 1999a,b). Al- 
though this had been checked from the beginning of this 
series of articles, it remains to be explicitly shown that 
the same approach, without new parameters, provides the 
correct cluster multiplicity and its observed evolution with 
redshift. 

Thus, the main goals of this article are to: 

- check the predictions for the properties of clusters of 
an analytic method developed in previous studies which 
can describe in a consistent way the mass functions of 
various objects (clusters, galaxies, Lya clouds, etc.) while 
making the connection with other properties of the density 
field (correlation functions, counts-in-cells, etc.). 

- evaluate by means of the same method the distortion 
of the CMB anisotropics induced by these X-ray clusters 
through the Sunyaev-Zel'dovich effect. 

- introduce a simple model which can reproduce the ob- 
served temperature - X-ray luminosity relation of groups 
and clusters. 

- use the global scope of our description to compare 
the X-ray emission provided by clusters, groups, galaxies 
and quasars. In particular, we show that this requires us 
to go beyond the standard PS prescription in order to deal 
with various classes of objects which may not be defined 
by the usual constant density threshold lS.c{z) ^ 177. 

In Sect|2| we present the expressions of the multiplic- 
ity functions and we compare the PS approach with the 
scaling model, focussing on the cluster mass function. The 
temperature function and its evolution with redshift are 
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discussed in Sect.^ while Sect.^is devoted to the implica- 
tions of our model for the characteristics of the Sunyaev- 
Zel'dovich effect. Next, in Sect.|| we present a model for 
the temperature - X-ray luminosity relation and we derive 
the cluster X-ray luminosity function. Finally, in Sect we 
describe the galaxy and quasar luminosity functions pre- 
dicted by our model, which allows us to draw a complete 
picture of the X-ray emission from all structures. 

2. Multiplicity functions 

2.1. Formalism 

We define clusters as halos with a mean density contrast 
equal to Ac, where Ac is the density contrast at the time 
of virialization, at the redshift z we consider, given by the 
spherical collapse model. This means that clusters just 
virialize at the redshift at which we see them. For a crit- 
ical density universe this density threshold is a constant: 
Ac ~ 177. Numerical simulations show that this value of 
the density contrast separates reasonably well the virial- 
ized halos from the surrounding material still falling onto 
the overdensity (Cole & Lacey 1996). This justifies this 
traditional definition of clusters. 

Following the method outlined in VS we shall use two 
prescriptions to get the comoving mass function of these 
halos. Firstly, we recall the scaling mode? developed in pre- 
vious papers, which we apply here to clusters of galaxies. 
Secondly, we also consider the usual PS approximation 
for the sake of comparison. Undoubtedly, we have in mind 
that the former supercedes the latter in the sense that it 
takes full benefit of the hierarchical clustering picture. The 
main advantages of this scaling approach are i) to make 
the link between the mass functions and the counts- in-cells 
statistics and ii) to provide a very powerful tool which can 
describe many different mass functions (i.e. defined by var- 
ious density thresholds) as well as other properties of the 
non-linear density field. 

The scaling model assumes that the many-body corre- 
lation functions follow specific scaling laws obtained from 
the stable-clustering ansatz (see VS for details). Then, we 
attach to each object a parameter x defined by: 

x{M, z) = _ , (1) 

e[i?(M,z),z] 

where 

m-J^^-^^Urur,) with V^jnR' 

is the average of the two-body correlation function 
'^2(ri, r2) over a spherical cell of radius R and provides the 
measure of the typical density fluctuations in such a cell. 
Thus large x correspond to deep, and small x to shallow 
potential wells. Then, we write the multiplicity function 
of these objects, defined by the density threshold Ac, as 
(see VS): 

V{M)— = ^ x^H(x) — , 



where Pq is the mean density of the universe at z = 0. The 
scaling function H(x) only depends on the initial spec- 
trum of the density fluctuations and must be obtained 
from numerical simulations. However, from theoretical ar- 
guments (see VS, Bernardeau & Schaeffer 1992 and Balian 
& Schaeffer 1989) it is expected to follow the asymptotic 
behaviour: 



X < 1 ; H{x) cx x"-^ 

with to ~ 0.5, ijJs ^ —3/2, a;» 
it must satisfy: 

/•OO 

/ X H{x) dx — 1. 
Jo 



a; > 1 : H{x) cx x'^ 



10 to 20 and by definition 



(3) 



This formulation is directly linked to the statistics of the 
counts-in-cells which involve a scaling function h{x)., re- 
lated to the probability P(A, i?) to have a density contrast 
A in a spherical cell of flxed radius _R, that scales as: 



(l + A)P(A,i?)dA = a;^ /i(a;)— , a;(A 



1 + A 

l{R,z\ 



(4) 



This function h(x) is very close to H{x), see VS and 
Valageas et al. (2000a) for details. The relevance of these 
scaling functions has been checked using numerical sim- 
ulations for CDM initial conditions by Bouchet et al. 
(1991) and, more systematically using power-law initial 
spectra in Colombi et al. (1997), Munshi et al. (1999) and 
Valageas et al. (20G0a). Here we note that in principle this 
prescription only applies to the highly non-linear regime 
(C ^ 100) while clusters correspond to mildly non-linear 
scales (^ '--^ 10 — 100). Hence the cluster mass functions 
we obtain may show an accuracy of 10%. Note however 
that the scaling function we shall use, taken from numeri- 
cal simulations by Bouchet et al. (1991), was measured in 
this range of ^, and was indeed found to describe cluster- 
ing at the above accuracy. 

On the other hand, the usual PS approximation gives 



7j{M) 



dM 



2 Po ^ 
TT M a 



diner 



dlnM 



exp 



2a2 



dM 



(5) 



where 6c{z) is the present linear density contrast for halos 
which collapsed at redshift z according to the spherical 
collapse model. As usual, a{M) is the rms density fluc- 
tuation extrapolated by linear theory at z = at scale 
M. Here we multiplied the mass function (||) by the usual 
empirical factor of 2 so that all the mass is contained in 
such overdensities. 

We can see that the scaling mass function (|^) predicts 
more numerous very massive halos but fewer small objects 
than the PS mass function (|^). This difference can be di- 
rectly seen through the scahng function H{x). Indeed, as 
shown in VS the PS formulation can be translated into the 
scaling approach in the highly non-linear regime where it 
leads to a specific scaling function: 



(2) hps{x) 



2 5 

TT 



6a 



n 5+2 
-X 6 



■ exp 



_^(5+n)/3/(2a2)j (6) 
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where n is the slope of the power-spectrum and a ~ 1. 
A comparison with the function h(x) directly measured 
in numerical simulations from counts-in-cells (see VS and 
Valageas et al. 2000a) shows that the large-mass cutoff 
of the PS mass function is too sharp while its peak (at 
masses where cr(M*) = 1) is too high. This latter re- 
gion corresponds to low mass clusters in Figj^ and to low 
temperature objects in Fig.^ and Fig.^ below. Note how- 
ever that at very low masses, corresponding to very small 
galaxies, the PS mass functions give fewer objects than 
the scaling prescription as seen in Valageas & Schaeffer 
(1999). 

Since in this article we wish to obtain in a consistent 
fashion the X-ray emission from all discrete sources: clus- 
ters, quasars and galaxies, we also need the multiplicity 
function of galaxies. As described in Valageas & Schaeffer 
(1999) we define galaxies by two constraints to be satis- 
fied simultaneously: the usual density threshold Ac{z) (as 
for clusters) and a cooling condition which ensures that 
the gas is able to dissipate its energy and form stars. This 
implies that galactic halos are defined by a density con- 
trast Agai{x, z) which depends on the mass of the object. 
Then, we can still use the expressions (|l|) and ^ to ob- 
tain the galaxy mass function, where the density contrast 
is now set to Agai{x,z) and depends on x (see VS and 
Valageas et al. 2000a). Moreover, from the identification 
(^ we can also consider an "extended PS" prescription to 
count these halos. We shall use this below in Sect.||when 
we deal with galaxies and quasars. 

For the numerical calculations we shall consider two 
cosmologies. First, we study a critical density universe 
(SCDM) with = 0.04, Hq = 60 km/s, erg = 0.5 and a 
CDM power-spectrum (Davis et al. 1985). Next, we con- 
sider an open CDM universe (OCDM) with = 0.3, 
Qa = 0, Ob = 0.03, Ho = 60 km/s and erg = 0.77. These 
values are those we used in previous articles where we con- 
sidered the luminosity functions of galaxies (Valageas & 
Schaeffer 1999), Lyman-a absorbers (Valageas et al. 1999) 
and reionization by stars and quasars (Valageas & Silk 
1999a). Thus, this present study of clusters and groups of 
galaxies completes our description of structure formation 
in the universe, so that we obtain a unified model which 
can describe in a consistent fashion all these objects, from 
small low-density Lyman-a absorbers up to massive clus- 
ters. 

Recent CMB observations (e.g., de Bernardis et al. 
2000) suggest that the universe is fiat. However, for our 
purposes, a ACDM universe with Sim = 0.3 and SIa = 0.7 
should be very close to our OCDM scenario, with a slight 
modification of the parameters of our astrophysical model 
for galaxies and clusters. In particular, as seen in Peacock 
& Dodds (1996) the functional relation between the linear 
and non-linear power-spectra is the same for the OCDM 
and ACDM (i.e. one simply needs to take into account 
the variation of the linear growth factor). Moreover, the 
scaling model which leads to the multiplicity function (H) 



has been seen to agree with the statistics of the projected 
density k along the line of sight obtained from numerical 
simulations for all three SCDM, OCDM and ACDM cos- 
mologies, as shown for instance in Valageas (2000). These 
results strongly suggest that our approach could also be 
used for the ACDM case. (In a similar fashion, the accu- 
racy of the PS prescription is similar for all three cosmolo- 
gies). 

2.2. Evolution with redshift of the cluster mass function 




S 1.5 1 0.5 O 0.5 

log(M/l O'f M^) 



Fig. 1. Upper panel: the cluster comoving cumulative mass 
function at redshifts z = 0.05, z = 0.33 and z = 1, for 
a critical density universe with a CDM power-spectrum. 
Higher redshifts lead to fewer massive clusters. The solid 
lines correspond to the scaling prescription while the 
dashed lines represent the PS formulation. Lower panel: 
same curves for an open universe with Sl,„ = 0.3. 

We display in Fig.|l| the cluster comoving cumulative 
mass function we obtain for the two cosmological scenarios 
we consider in this article, at redshifts z — 0.05, z = 0.33 
and z = 1. First, we can check that we recover the trend 
described in the previous section: the scaling approach 
predicts fewer low mass halos and more very massive ob- 
jects than the PS prescription. In a similar fashion, Gov- 
ernato et al. (1999) find that the PS mass function overes- 
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timates the number of small clusters (M < 2 for 
(jg = 0.7) and underestimates the number of massive halos 
(see also Gross et al. 1998 for a similar trend). However, 
the discrepancy they measure is smaller than ours. This 
may be due to the fact that the power-spectrum we use 
(from Davis et al. 1985) is slightly different from theirs 
and to the use of our scaling model in a range of ^ slightly 
beyond its range of validity: indeed very high masses cor- 
respond to large scales which are getting close to the linear 
regime where the scaling model is not valid. Nevertheless, 
the agreement of our predictions (already described in a 
more general context in VS), with the behaviour observed 
in the numerical simulations, is quite encouraging. 

We can check that the redshift evolution of the mass 
function is very sensitive to the cosmological parameter 
rim (e.g., Oukbir & Blanchard 1992; Eke et al. 1996). In- 
deed, the number of clusters declines faster with z for the 
critical density universe than for the open model. This is 
simply due to the fact that in the latter case structures 
have nearly stopped growing since the redshift where 
became appreciably smaller than unity, whence a very 
weak evolution with redshift, while in the former case 
structures keep building under the action of gravity at 
all epochs. 



where dV/dfldz is the comoving volume element per unit 
steradian and unit redshift. The normalization of the PS 
prediction is lower than for the scaling approach as ex- 
plained above (we count halos which are in the tail of 
the mass function). The main result of this figure is to 
emphasize the difference between the two cosmologies of 
the redshift evolution. Thus, for ri,„ = 1 the number of 
such clusters reaches a maximum at z 0.3, while for 
f2,„ = 0.3 the peak corresponds to z ~ 0.45 and the evo- 
lution is slower. 

2.3. Counts along the line of sight 



o 




o 




-a -1.5 -1 -0-5 O 0-5 

log(M/l Qio Mg) 

Fig. 3. The mean number of clusters Afis{> M) of mass 
larger than M which intersect a line of sight between z = 
and z — 2 (the contribution of higher redshifts is negligi- 
ble). The solid lines correspond to the scaling prescription 
and the dashed lines to the PS approach, for both SCDM 
and OCDM scenarios. In both cases the number of very 
massive clusters is larger for the open universe. 



O 0.^^ 0-4 0.6 O.B 1 

z 

Fig. 2. The redshift distribution per square degree of 
clusters more massive than 5 x 10"'^^ M©, for fim = 1 
and l^m = 0.3. The graph shows the number of clusters 
dM /dz{> 5 10^'* Mq) per unit redshift interval and square 
degree. The solid lines correspond to the non-linear pre- 
scription, while the dashed lines represent the PS formula- 
tion. The open universe corresponds to the slowest redshift 
evolution. 

The dependence on of the redshift evolution of the 
cluster mass function is even more apparent in Fig.^ which 
shows the redshift distribution of clusters more massive 
than 5 x 10^'* Mq, per square degree, for both SCDM and 
OCDM scenarios: 

dU , 14 , f n \ dV , ,dM , , 

-O5X10-M,). (-) viM)-il) 



From the cluster multiplicity function we can derive 
the mean number of clusters Afis{> M, < z) of mass larger 
than M which intersect a line of sight between z = and 
a given redshift z: 

dt dM 

Msi>M,<z)= Al + zfdz ^^{M)nR\ 
Jo dz Jm M 

The result is displayed in FigJ|. Of course, we recover the 
difference between the scaling model and the PS approach 
we described above for the mass functions. The open uni- 
verse gives higher counts for very massive halos because of 
the slower decline with redshift of the mass function. Note 
however that for the scaling prescription the difference be- 
tween the two cosmologies is quite small. Of course this 
would change with another choice for We see that the 
mean number of clusters on a line-of-sight is quite small 
since we typically have Nis < 1. 
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3. Cluster temperature function 

Although the study of cluster mass functions is convenient 
from a theoretical point of view, for observational purposes 
it is more interesting to consider temperature functions. 
For this we need the temperature which is associated with 
halos of a given mass. 

3.1. Characteristic temperature of the halos 

The Jeans equation for the velocity dispersion of the 
dark matter yields: 



d 
dr 



(P'^l) = - P 



gM(< r) 



(8) 



In the case of an isothermal density profile p{r) oi r ^ this 
leads to: 



''viR)^^ and kT- 



Q^mpM 
2R ' 



(9) 



where we defined the virial temperature T of the halo by 
kT = firupa^ and R is the virial radius of the cluster. Here 
lirrip is the mean molecular weight of the gas and mp is the 
proton mass. For a different density profile we would still 
obtain (^) for the mean temperature, with a multiplicative 
factor of order unity. We assume that the mass of baryons 
Mb is proportional to the mass of the dark matter halo 
M: 



Mb = ^ M, 



(10) 



where fib is the present ratio of the baryon density to the 
critical density. With these parameters, using the fact that 
halos are defined by the density contrast Ac(z), we have: 



(11) 



where R is the virial radius of the cluster, and eq. 
writes: 



T = To Aff/' Ac(z)'/' (1 + z) 
with 

M 



(12) 



Ml 



1015 Mp 



^ and To = 1.2 /i^/s ]^^y 



This is consistent with numerical simulations which re- 
cover this scaling law, with a similar normalization. In 
these units, Navarro et al. (1995) find To = 1.4 ft-^/s 
keV while Evrard et al. (1996) get Tq = 1.2 /i^/s 

3. 2. Evolution with redshift of the temperature function 

Using ri{T)dT/T = r]{M)dM/M we can get the cluster 
temperature function from ( p^ ) and (^ or (0). Its evolu- 
tion with redshift is shown in Fig.^ and Fig^ Note that 
the temperature T we consider in this section is the virial 
temperature. Indeed, for the hot clusters (T <; 1 keV) wc 




Fig. 4. Upper panel: the comoving cluster X-ray temper- 
ature function r]{T) dT/T for the SCDM scenario. The 
solid lines correspond to the scaling prescription while the 
dashed lines represent the PS formulation. We display the 
redshifts z = 0.05, z — 0.33 and z = 1 {a, larger z corre- 
sponds to fewer bright clusters). The data points are ob- 
servational results at z = from Henry & Arnaud (1991) 
(disks) and Edge et al. (1990) (squares). Lower panel: the 
cumulative cluster temperature function at the same red- 
shifts. The open and filled squares are observations from 
Henry (1997) at 2 = 0.05 and z = 0.33 respectively. 



study here it is also the temperature of the gas which is 
heated by shocks during the gravitational collapse. In con- 
trast, in cool groups {T <, 1 keV) the gas is also influenced 
by a possible preheating of the IGM (e.g., Valageas & Silk 
1999b) which leads to a smoother baryonic density profile 
and a larger gas temperature. This is discussed in Sect.|^. 

First, we can check that the difference between the 
scaling prescription and the PS prediction is similar to 
the trend we obtained for the mass functions in Sect. ^.2[ . 
In particular, we can note that Governato et al. (1999) 
found that for a standard SCDM model {flni = 1) normal- 
ized to (Tg = 0-5 (which is our case) the PS prescription 
underestimates the number of hot clusters kT > 7 keV 
by almost a factor 10 at z = 1, while it overestimates the 
number of small halos. These authors also found that the 
deficiency of massive halos predicted by the PS approach 
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Fig. 5. The comoving differential and cumulative temper- 
ature functions as in Figj^, for the case of an open universe 
with rim = 0.3. 



gets more severe for smaller as (hence also at higher red- 
shifts). In our view (see the discussion in VS), this is sim- 
ply due to the fact that in this case one looks at rarer 
objects, farther in the cutoff of the mass function, which 
increases the discrepancy between both theoretical mass 
functions which have different exponential tails. However, 
the discrepancy between the PS and the scaling predic- 
tions is again slightly larger than what these numerical 
results imply at the large mass end. 

Next, we can check that the redshift evolution of the 
temperature function is faster for the critical density uni- 
verse than for the open case, in agreement with Fig.|l| for 
the mass function (see also Oukbir & Blanchard 1997). 
Thus, the very small decline with z of the observed clus- 
ter temperature function (Henry 1997) from z = 0.05 to 
z = 0.33 favors the open case (or more generally a low- 
density universe). However, the redshift evolution of the 
cluster temperature function we obtain in the SCDM case, 
is not much faster than the observed decline and our study 
shows it cannot be ruled out. Note that we could obtain a 
better agreement with the data for ri,„ = 1 by choosing a 
slightly lower cg. However, the normalization of the power- 
spectrum we use is constrained by our previous studies of 
galaxies, quasars and Lyman-a clouds since we want to 




Fig. 6. The redshift distribution per square degree of clus- 
ters hotter than 5 keV, for Q.^ = 1 and = 0.3, f^A = 0, 
with a CDM power-spectrum. The graph shows the num- 
ber of clusters d/V /dz{> 5 keV) per unit redshift interval 
and square degree. The solid lines correspond to the non- 
linear prescription while the dashed lines represent the PS 
formulation. The open universe corresponds to the slowest 
redshift evolution. 



build a unified consistent model. Hence we must choose 
a value which provides good results for all these objects. 
Moreover, as explained in Sect.|^ massive clusters corre- 
spond to mildly non-linear scales close to the theoreti- 
cal limit of validity of the scaling model, so that we may 
slightly overestimate the number of very massive clusters. 
Thus the value trg = 0.5 used for = 1 (which is also 
the result obtained by Governato et al. 1999) seems satis- 
factory. 

The redshift evolution of the temperature function is 
slower than the change in the mass function which was 
presented in Figj^. This is due to the temperature-mass re- 
lation (^ which imphes that T cx M^/^ A^izf/^ (1 + z). 
Thus, the temperature which corresponds to a given mass 
increases with z, which enhances the redshift evolution 
of the mass function as compared with the temperature 
function. 

The redshift distribution of clusters hotter than 5 keV, 
per square degree, is shown in Fig|6[ Of course, our results 
are similar to Figj2[ 



4. Sunyaev-Zel'dovich effect 

An indirect method to get observational constraints on the 
cluster temperature function is to measure the Sunyaev- 
Zel'dovich (SZ) effect (Sunyaev & Zel'dovich 1972) due 
to the hot gas in the intra-cluster medium. As compared 
to the X-ray luminosity function which we discuss below 
in Sect.||, the SZ effect presents two strong advantages: 
it does not depend on the detailed density profile of the 
gas distribution within clusters and it is more sensitive 
to high-redshift objects. Hence in this section we describe 
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our predictions for this indirect measure of the cluster 
multiphcity function. 

First, we recall that the variation of the CMB bright- 
ness at the frequency v along the line-of-sight due to the 
SZ effect can be written as (e.g., Barbosa et al. 1996): 



(13) 



where jv describes the spectral form of the distortion, in- 
dependent of the cluster, and y is the Compton parameter 
given by an integration along the line-of-sight through the 
cluster: 



V 



UeCTT- 



(14) 



Here, Tg is the temperature of the electrons in the intra- 
cluster gas (which we approximate by the virial temper- 
ature T), me the electron mass, Ue the electron number 
density and ar = 6.65 10~^^ cm^ the Thompson cross sec- 
tion. Defining the dimensionless frequency u = hpv/kTo = 
Ao/A where hp is Planck constant and Aq = 5.28 mm for 
To = 2.726 K, To being the present temperature of the 
CMB, one can write: 



Mu) = 2 



(kTo) 



(hpc)^ (e" - 1)2 Vtanh(u/2) 



-4 



(15) 



The flux S^{u) of the cluster, in mJy = 
10~26 erg s~^cm~2Hz~"'^, is simply the integral of ii, 
over the solid angle subtended by the cluster: 



S^{u) = j^{u) rd{z) ^ 



kT 

UeCFT ^dV 

meC^ 



(16) 



where: 



rd{z) 



2c 



Honiii + zy 



(fin, - 2) + ^n,Z - 1 



is the angular distance of the cluster, located at redshift z. 
Hence, the total flux observed from an unresolved cluster 
depends only on the mass of gas at the temperature Tg, 
and not on the density profile. Moreover, we can see that 
y only depends on the physical properties of the cluster, 
and not on its redshift. Hence it is very sensitive to the 
cluster populations at high redshifts, which contribute in 
the same manner as close clusters, which means it is a 
useful tool to study the distant universe. Finally, in the 
long wavelength regime {u — > 0) the negative fluctuation 
of the CMB spectrum is simply given by: 



AT 

Y 



-7^ = -22/. 



(17) 



The total mean Compton parameter (y) , averaged over 
all lines of sight, which only depends on the temperature 
distribution of the gas, can be obtained from ( p^ : 



(y) 



cdt / <tt 



kTg drie ^rj. 

nieC^ dTq ^ 



(18) 



<y> 












<y> 












OCDM 






n,^= 0.3 





Fig. 7. Upper panel: the Compton parameter y, and its 
fluctuations (Sy"^)^/^, on a line-of-sight from z = up to 
the redshift z, for the case f^m — 1- The solid lines corre- 
spond to the scaling formulation and the dashed lines to 
the PS prescription. In both cases the fluctuation (5y2^i/2 
is the curve which shows the steepest rise at z = 0. Lower 
panel: same curves for an open universe with fl^ — 0.3. 

where Ue is the mean electron number density. Whence: 

f°° dt 
(y) ^ / dz—caril + z)^ 
Jo dz 



Mi 



dM , , f7b M kT„ 

M Sim MlTOp rUeC^ 



(19) 



The factor (1-1- 2)'^ transforms the comoving mass function 
r]{M) we have used so far, into the mass function in proper 
coordinates. The cutoff at Mi{z) (defined by tcooi < tn) is 
due to the fact that in small halos the gas cools in less than 
one Hubble time tn- The mass Mi{z) we get is typically 
of the order of lO^^ m^^^ 

An alternative way to get (y) is to divide the line-of- 
sight into small length elements A/^, so that: 



{y)^{Y.y^) and {y')^dY. 



(20) 



where iji is the Compton distortion due to the line element 
A/i in one random realization. In this way, we recover the 



Valageas & Schaeffer: Multiplicity Functions and X-ray Emission of Clusters and Groups versus Galaxies and Quasars 9 




III 


1 1 1 1 1 1 1 

OCDM 


III 




= 0.3 












1 1-5 2 2,5 3 

log(S„) [ mJy ] 

Fig. 8. Upper panel: the SZ source counts at A = 0.75 mm 
for the SCDM case. The sohd hne corresponds to the scal- 
ing formulation and the dashed line to the PS prescription. 
Lower panel: same curves for an open universe. 



previous formula for (y) and we get the fluctuations 
of y: 



(V>-((y-(y)n 

with: 



(21) 



f°° dt 
(5y^) = J dz—caxil + z) 



R? ^iTTlp TTleC^ 

+aT£.cc{Rcuz)Rci{z){l + zf 



(22) 



Mi 



dM M kTg 

M r2,n iiirUp m^c^ 



The first term is due to the Poisson fluctuations of the 
number and the mass of clusters in each line element 
while the second term arises from the correlations $,cc{f', z) 
among clusters, and Rci{z) is the typical radius of clusters. 
The first term scales as the number density Ud of clusters 
while the second term scales as nci{£,cc{R)nciR^ci) ■ Thus, 



the Poisson fluctuations are greater since the number den- 
sity of clusters is small {£,cc{R)nciR% ^ 1). We show in 
Fig.0 the mean Compton parameter (y) and its fluctu- 
ations {Sy"^)^^^ on a line-of-sight from z = up to the 
redshift z. The Compton parameter (y) is larger for the 
open universe because the cluster mass function declines 
more slowly at higher z and the line-element is slightly 
larger. Since the number of clusters on the line-of-sight is 
rather small (see Fig.^) the fluctuations of y are of the 
same magnitude as the mean (y). Of course (5y^)^/^ is 
much larger than (y) for z ~ since the number of clus- 
ters tends to in this limit. In particular, at z ~ we 
have: 



(y) cx [(1 + z)3/2 - l] and {SyY^"^ cx (y)!/^. 



(23) 



Note that the COBE/FIRAS upper limit is (y) < 1.5 x 
10^^ (Fixsen et al. 1996). Thus, in both cosmologies the 
Sunyaev-Zel'dovich effect as measured by y is still one 
order of magnitude below the upper limit provided by 
present observations. Such a distorsion, however, will be 
within the reach of the MAP and PLANCK projects. 




Fig. 9. The redshift distribution of SZ source counts for 
the SCDM and OCDM scenarios at A = 0.75 mm. The 
solid lines are the scaling formulation and the dashed 
lines the PS prescription. The low-density cosmology cor- 
responds to the slowest decline at large z of the SZ counts. 

From the cluster mass function we can derive the flux 
density distribution of clusters i]{Si^)dSv / S^,- Hence, the 
number of SZ sources per unit solid angle on the sky with 
a total monochromatic flux greater than Si, (we assume 
that the objects are unresolved) is 



— > 5^) = / dz—— / r](M,z)—— 
dn^ y dfldz Jm' ' M 



(24) 



The cutoff M[{S,j^ z) corresponds to the threshold S^- The 
SZ source counts at A = 0.75 mm are shown in Fig.|| as 
a function a 5*1, for both cosmologies. As for the Comp- 
ton parameter y, the open universe leads to slightly higher 



10 Valageas & SchaefFer: Multiplicity Functions and X-ray Emission of Clusters and Groups versus Galaxies and Quasars 



counts because of tlie slower decline of the cluster multi- 
plicity function and of the larger volume element. 

The redshift distribution of these source counts is dis- 
played in FigJ|. Of course, as explained above, the open 
universe predicts a slower decline at large z of the source 
counts than for a SCDM scenario. 

5. Evolution of the X-ray luminosity function 

Although observations can provide an estimate of the clus- 
ter temperature function they can more easily give the 
X-ray luminosity function, since it is easier to measure 
the luminosity of a distant cluster than its temperature. 
Hence we describe in this section a model to obtain the 
X-ray luminosity of clusters of galaxies. Moreover, the 
temperature-luminosity relation also contains some inter- 
esting information on a possible reheating of the IGM. 

5.1. The temperature-luminosity relation 

5.1.1. Breaking the simple scale-invariance 

The bolomctric X-ray luminosity L^oi of a cluster of vol- 
ume V is 



^bol 



nlKb{Tg)dV 



(25) 



where e is a constant of order unity, is the electron 
number density and Ki)[Tg) is the emissivity function (in 
erg cm'^ s~^) for a gas at the temperature Tg. Thus, con- 
trary to the Sunyaev-Zel'dovich effect, see (|l6|), the X- 
ray luminosity strongly depends on the density profile of 
the hot gas within the cluster. Using (|l2|) and the fact 
that Ab(T) c>c \/T for bremsstrahlung one expects (Kaiser 
1986): 



(26) 



Note that this power-law behaviour is related to the scale- 
invariance of clusters: more massive objects are identical 
to smaller ones after a simple rescaling. However, obser- 
vations show a much steeper slope Lbo; T^-^^ (Arnaud 
& Evrard 1999) and suggest a bend in the temperature- 
luminosity relation. Hence some physics is missing in the 
derivation leading to (p6|). In particular, one needs to 
break the scaling laws which led to (|2^) through the intro- 
duction of additional dimensional quantities. It has been 
suggested in the litterature (e.g., Evrard & Henry 1991; 
Cavaliere et al. 1997; Ponman et al. 1999; Lloyd-Davies et 
al. 2000) that one needs to take into account the reheat- 
ing of the IGM, prior to cluster formation, which raises 
its entropy and can modify the gas dynamics in cool clus- 
ters. Indeed, if this entropy "floor" is sufficiently large, the 
gas is heated by the adiabatic compression up to a tem- 
perature which can be as large as the virial temperature 
of the halo. In this case the density profile of the gas is 
much smoother than the distribution of the dark matter. 



which diminishes the luminosity of these cool clusters and 
modifies the relation Tg — Lf,oi. Such a reheating of the 
IGM by quasars or supernovae was studied in Valageas & 
Silk (1999b) where it was shown that the energy provided 
by quasars may be sufficient to reheat the IGM. The en- 
ergy delivered by supernovae was found by the latter au- 
thors to be rather small, but a more thorough estimate 
is under way. In this article, we simply assume that such 
processes have reheated the gas to a characteristic tem- 
perature Tad ~ 0.4 keV. This breaks the simple scaling 
laws and can lead to a non-trivial temperature-luminosity 
relation. 

5.1.2. Density profiles of the gas and of the dark matter 

In order to obtain the distribution of the gas, we need the 
density profile of the dark matter, which is gravitationally 
dominant. Here we assume that the halos obey the density 
profile obtained in high-resolution simulations by Moore 
et al. (1999b): 



p{r) = 



(r/r,)i-5 [1 + (r/r,)i-5] 



with 



R 

rs = — 
c 



(27) 



where c is the concentration parameter. We use c — A 
as in Moore et al. (1999b), although this parameter may 
slightly depend on the mass of the cluster. Moreover, we 
assume hydrostatic equilibrium for the gas. We consider 
two models for the temperature profile Tg{r) of the gas. 
First, we write: 



Tg{r)^T,{r)+Tad 



(1 + A, 



(28) 



where 7s — 5/3. The first term on the right hand side in 
( p8| ) describes non-adiabatic gravitational heating through 
shocks during the formation of the cluster. The second 
term takes into account the reheating of the gas, before 
the formation of the cluster, and its subsequent heating 
through the adiabatic compression during the infall. Note 
that the relation ( p8| ) assumes that the thermal conduction 
is small (e.g., Sarazin 1988) so that the gas is not exactly 
isothermal. In a fashion similar to (O) we write: 



Ts{r) - 



Q^mpM(< r) 
2kr 



ln(l + (r/r.)i-^) R 

ln(l-hci-5) r ^ ' 



which also measures the depth of the potential well. This 
also satisfies the Jeans eq.(||) within a factor 2. Then, the 
distribution of the gas is given by the condition of hydro- 
static equilibrium: 



dP_ 

dr 



Pg 



gM{< r) 



with P 



/imp 



(30) 



However, in practice we use a simplified model. From (30) 
one can check that at large radii where Tg ~ Ts{r) the 
gas follows the dark matter while at small radii below 
Rcore where both terms in (E8|) are equal, the gas density 
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saturates. Indeed, when Tg ^ Ts the gas no longer falls 
into the potential well and a core with a nearly constant 
baryon density appears. This model is similar to those 
used in Cavaliere et al. (1997,1999) or Wu et al. (1999) 
except that these authors use the dark matter density pro- 
files given by Navarro et al. (1996) and they do not take 
into account the first term on the right hand side in ( psf ) 
(they use a polytropic equation of state). However, in all 
cases the behaviour of the gas distribution is the same: for 
smaller halos the preheating becomes more important and 
the gas density profile gets flatter which decreases the X- 
ray luminosity. Here, in order to simplify the calculations 
we approximate the gas distribution by the same profile 
( p7| ) as for the dark matter in the outer parts r > Rcore 
and by a constant value in the core r < Rcore'- 



r> Rc 



r < Rc 



p{r) 



(31) 



which ensures that the total mass is conserved. 

Second, we consider an alternative model where the 
temperature of the gas is given by: 



(32) 



As in ( Pq ) two processes govern the temperature of the gas: 
gravitational heating (first term) and a second unspecified 
source of energy (e.g., SNe or QS Os) which breaks the scal- 

(second term). However, 



5.1.1 



ing law described in Sect, 
in contrast to (Pa), here we take this additional term to 



be constant. Hence this model corresponds to a uniform 
"post-heating" of the gas which occurs after the forma- 
tion of the cluster, for instance through the SNe which 
may eject some energy into the intra-cluster medium dur- 
ing and after the infall of the gas into the potential well 
of the cluster. From another point of view, this can also 
be seen as a specific case of the "pre-heating" model ( p8|) 
with a different equation of state: 7s ~ 1. Then, this al- 
lows us to estimate the sensitivity of our model (pq) to 
the assumed value for 7s. Next, from the relation ( |32| ) we 
proceed exactly as from (^) to obtain the density distri- 
bution of the gas, described by a core radius Rcore as in 
(^), which is now defined by Ts{Rcore) — Tad- 

5.1.3. Cooling radius 

In addition to the core described above which appears for 
the gas distribution, cooling may affect the distribution of 
hot X-ray emitting gas. Indeed, in the inner parts of the 
cluster the density can be large enough to lead to a small 
cooling time so that a cooling flow develops. Then, some 
of the gas forms a cold component which does not emit in 
X-ray any longer. Thus, we define the cooling radius Rcooi 




log(T ) 



Fig. 10. The relation temperature - bolometric luminosity 
for clusters in the case 7s = 5/3 for the model ( [2^ ) with 
Tad = 0.5 keV (resp. Tad = 0.4 keV) for the SCDM (resp. 
OCDM) scenario. The solid, dashed and dot-dashed lines 
correspond to z = 0, z = 0.33 and z — 1 respectively. The 
data points are from Mushotzky & Scharf (1997) (0.14 < 
z < 0.55) for clusters and from Ponman et al. (1996) {z < 
0.05) for groups. The upper panel corresponds to a critical 
density universe (SCDM) and the lower panel to an open 
universe (OCDM). 



as the shell where the cooling time tc 
the Hubble time t^iz): 



becomes equal to 



tr 



tn with tc 



2npgAc{Tg) ' 



(33) 



where Ac{Tg) is the cooling function. At large radii r > 
Rcooi the local cooling time is larger than the Hubble time. 
In this case, the gas distribution and the temperature have 
not had time to evolve much and the X-ray emissivity is 
proportional to Ab{Tg), see (p5|). On the other hand, 
within the cooling radius Rcooi the gas has had time to 
cool and to form dense cold clouds. However, we consider 
that some of the gas is still hot and emits in X-ray as cool- 
ing does not proceed in a uniform fashion (Nulsen 1986; 
Teyssier et al. 1997; Waxman & Miralda-Escude 1995). 
Indeed, the cooling instability leads to a wide range of gas 
temperatures and densities as overdense regions cool faster 
and contract (because of the pressure of the surrounding 
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Fig. 11. The relation temperature - bolometric luminosity 
as in Fig.|l^ for the model (32) (i.e. 7^ = 1) with the 
same values for Tad- A larger value for Tad would improve 
the agreement with observations at low Tg for the OCDM 
scenario. 



gas) which increases even further their density contrast. 
Then, some of the gas is simply removed from the X-ray 
emitting component as these dense cold clouds decouple 
from the hot phase. On the other hand, the temperature 
of the hot gas remains close to Ts{r), introduced in (p9|), 
through adiabatic compression and possible gravitational 
heating, see Nulsen (1986, 1998) for detailed models of 
this multiphase ICM. The density of the hot component 
must be of order pcooi, defined by the condition tcooi = tn, 
since the density has not had time to decrease further yet 
(note that the system is not stationary). As this multi- 
phase medium is connected to the outer parts of the clus- 
ter r > Rcooi which also provide a reservoir of matter 
and energy, we assume that the hot low-density phase is 
spread all over the radius Rcooi- Note that the time-scale 
for hydrostatic equilibrium is ip ^ r/cs ~ tdyn for gas 
at the temperature Tg (where Cg is the sound speed while 
tdyn is the dynamical time) . Since the time-scale of the in- 
fall of the gas cannot be smaller than tdyn and is actually 
expected to be larger (the pressure and possible energy in- 
jection from SNe or gravitational interactions slow down 
the motion) there should be approximate pressure equi- 



librium. Hence at small radii r < Rcooi we write for the 
density of the X-ray emitting gas: 

r < Rcooi ■ Pgx{r) = Pcool = PgiRcool) < Pgir) (34) 

while at large radii r > Rcooi we have Pgx{r) — Pgir)- 
For cool clusters where the non-gravitational energy term 
Tad plays an important role, the distribution of the gas is 
flatter than the dark matter and it shows a core of constant 
density Pcore- Then, the cooling time at this core radius 
Rcore is still longer than the Hubble time so that cooling 
plays no role: the luminosity of the cluster is determined 
by the radius Rcore due to the distribution of the gas itself. 
On the other hand, for hot clusters the gas follows the 
dark matter density profile over a large range of radii and 
the cooling time gets smaller than the Hubble time in the 
"outer region" r > Rcore- In this case, the distribution 
of the X-ray emitting gas is characterized by the cooling 
radius Rcooi and the luminosity of the cluster is governed 
by the cooling criterium. Note that Rcooi mainly plays the 
role of a cutoff for the distribution of X-ray emitting gas: 
we would obtain similar results for a model where we set 
Pgxir) = for r < Rcooi, as most of the X-ray emission 
comes from the regions r ~ Rcooi (thus such a model would 
simply decrease the luminosity by a numerical factor ~ 2 



which can be absorbed in the normalization of tr 



for 



instance). 

In Valageas & Silk (1999b) we used for illustrative pur- 
poses an isothermal model for the dark matter and the gas 
distribution, with the relation (|3|) (interpreted as 7s = 1). 
Thus, we had: 



p{r) cx r ^ and Pg{r) oc p ^ with (3 = 



T 



T + Tad 



, (35) 



where T is the virial temperature of the cluster defined in 
(^ and the gas distribution was obtained from the hydro- 
static equilibrium condition (^). Here, the distribution of 
the gas does not show a core radius Rcore while the X-ray 
emitting gas is still characterized by a cooling radius Rcooi - 
However, it is easy to check that we recover a behaviour 
similar to the case discussed above. Indeed, for hot clus- 
ters with T > 3Tad the density profile of the gas is very 
steep (/3 > 3/4) so that the X-ray emissivity is dominated 
by the inner parts of the cluster - as can be seen from (25) 
and ( p5[ ) - hence by the cooling radius Rcooi - On the other 
hand, for cool clusters with T < 3Tad the density profile 
of the gas is rather flat (/3 < 3/4) which means that the 
luminosity of the cluster is governed by the outer parts 
r ~ i?. In this paper, for cool clusters the emissivity is 

dominated by the regions T ^ Rcore S-Ild Rcore 

R for 

low temperature clusters (T < Tad)- 

Thus, we see that the main characteristics of clusters 
do not strongly depend on the details of the models (den- 
sity proflle of the dark matter halo, index 7^ of the "equa- 
tion of state") as they are mainly sensitive to Tad- 
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5.1.4. Evolution of the temperature - X-ray luminosity 
relation 

For both models described above we obtain the X-ray lu- 
minosity from (p5|), where rie is given by the density of 
the X-ray emitting gas pgx ■ We determine the factor e by 
requiring (p5| ) to reproduce the observations for massive 
clusters {Tg > 1 keV). We get e = 2.2 (resp. e = 0.8) for 
a critical density universe (resp. an open universe). We 
show in Fig.|l^ and in Figjl^ the temperature-luminosity 
relation we obtain for both models ( p^ ) and (^^, with 
the above value of e and Tad = 0.5 keV (resp. 0.4 keV) 
for the SCDM (resp. OCDM) cosmology The gas tem- 
perature used in the figures is Tg = T + Tad- We can see 
that we get very similar results which agree reasonably 
well with observations. Of course, we could improve the 
agreement of the model ( |3^ ) (i.e. 7^ = 1) with observa- 
tions in the OCDM cosmology by using a slightly larger 
value for Tad- We also note that the redshift evolution we 
obtain is very small, which is consistent with observations 
(Mushotzky & Scharf 1997). This suggests that the "pre- 
heating" temperature Tad should not evolve too strongly 
with redshift for z < 1, which might favour supernovae as 
the source of energy (due to the sharp decline at low z 
of the quasar luminosity function, the characteristic tem- 
perature Tad obtained within the framework of a model 
where the preheating is due to QSOs, is expected to show 
a faster evolution with z, see Valageas & Silk 1999b). 

Thus, in order to recover a bend in the temperature - 
X-ray luminosity relation one mainly needs to introduce 
a new dimensional parameter, like Tad, as discussed in 
Sect. I . 1 . l| . Since in this article we are only interested in the 
total luminosity of clusters the simple models described in 
Sect. I .1.2| are sufficient for our purpose and in the follow- 
ing we use the model ( p8| ) (i.e. 7^ = 5/3). As shown by the 
comparison of Fig. and Fig.^our predictions for the lu- 
minosity function are not very sensitive to the details of 
our model. We would also obtain similar results with the 
isothermal model ( ^5| ) as in Valageas & Silk (1999b). On 
the other hand, this means that in order to discriminate 
between various models one needs precise observations of 
the density profiles of the gas and of the dark matter as 
well as a measure of the gas temperature. This will be 
provided by the XMM mission. However, present observa- 
tions seem to favour the model (|2^) with a "preheating" 
of the gas before the formation of clusters (Lloyd-Davies 
at al. 2000). 

5.2. Luminosity Junction 

From the temperature-luminosity relation described in the 
previous section and the temperature multiplicity func- 
tion obtained in Sect.^ we can derive the cluster X-ray 
luminosity function. Note that the results discussed in 
the following are largely independent of Sect. ^.l] since any 
temperature-luminosity relation which agrees with obser- 
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Fig. 12. The comoving cluster luminosity function in the 
rest-frame 0.5 — 2 keV band, at redshifts z — z — 0.33 
and z = 1 (smaller z corresponds to fewer faint ob- 
jects). The solid lines show the scaling formulation and 
the dashed lines the PS prescription. The data points at 
z = are from Ebeling et al. (1997) (filled squares) and 
Burns (1996) (crosses). The upper panel corresponds to 
the SCDM scenario and the lower panel to the open uni- 
verse. 

vations would give similar results. From the luminosity 
Lfjoi obtained in the previous section we write the lumi- 
nosity Lx in the frequency band vi — vi as: 

Lx = Uoi [e-^^'^'l^^^ - e-''^"'2/fe^s) . (36) 

Here hp is Planck constant and we neglected the varia- 
tion of the Gaunt factor. We compare our predictions with 
observations in Fig.^ for both cosmological scenarios, in 
the rest-frame frequency band 0.5 — 2 keV. First, we note 
that we recover the fact that the scaling model predicts 
more massive and bright clusters but fewer small and faint 
objects than the PS approach. Then, we see that for the 
critical density universe, the luminosity function we obtain 
predicts too many faint clusters. This could be cured by a 
small change of the initial power-spectrum. Of course, in 
a similar fashion one can also bring the PS prescription in 
agreement with observations. However, we prefer to keep 
this normalization of the power-spectrum in order to be 
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consistent with our previous articles about galaxies and 
Lyman-Q! clouds and with results from numerical simula- 
tions (Governato et al. 1999). On the other hand, for the 
open universe our predictions agree reasonably well with 
the data. We can see that for both cosmological scenar- 
ios the redshift evolution we get is very small, which is 
consistent with observations. 




Fig. 13. The number of clusters per square degree brighter 
than an X-ray flux limit Sx- The solid line corresponds 
to the scaling formulation and the dashed lines to the PS 
prescription. The low dot-dashed curve is a no-evolution 
model. The data points are from Jones et al. (1998) (tri- 
angles) and Rosati et al. (1998) (disks). The upper (resp. 
lower) panel corresponds to the SCDM (resp. OCDM) sce- 



We display in Fig.|l3| the integrated cluster surface den- 
sity above an X-ray flux threshold Sx in the frequency 
band 0.5 - 2 keV: 



dJ^^ r, ^ f . dy r ^dM 



(37) 



The cutoff M[{Sxi z) corresponds to the X-ray flux Sx - It 
is obtained from the temperature-luminosity relation and 
the flux-luminosity relation: 



where the distance rium{z) is the luminosity distance up 
to redshift z. Here the luminosity Lx is obtained as in (36) 
with i'i{z) = + z) to take into account the redshift of 
the observed frequency band 0.5 — 2 keV. Our results agree 
with observations, although our normalization is slightly 
too high for the critical density universe. In both panels, 
the dot-dashed curve corresponds to a no-evolution model 
where the comoving cluster luminosity function does not 
vary with z and remains equal to its value at 2 = 0, shown 
in Fig.|l^. Note that evolution effects are not very large. 
The scaling and PS approaches give very close results, al- 
though we can still recognize that the PS approximation 
predicts fewer massive objects but more small halos. In- 
deed, this difference is somewhat smeared out in FigJT^ 
because the X-ray sources seen with a given flux Sx cor- 
respond to a large variety of objects located at different 
redshifts. 



6. Galaxies and quasars versus groups and 
clusters 

6.1. Galaxy X-ray luminosity Junction 

In addition to clusters, galaxies may also emit in X-rays 
when they form. Indeed, in order to make stars and build a 
galaxy the gas embedded within a dark matter halo must 
cool and fall into the gravitational potential well. During 
this process, the gas can radiate some energy in the X-ray 
band by bremsstrahlung, especially for the most massive 
galaxies with a large virial temperature T ^ 10^ K ^ 1 
keV. As a consequence, some of the X-ray sources one 
could observe on the sky may be high-redshift newly-born 
galaxies. Note that there will be obvious observational dif- 
ficulties to distinguish X-ray emitting galaxies from small 
groups of similar mass, containing a few smaller galaxies. 

To derive this galaxy X-ray luminosity function that 
is to be compared to the one for groups and rich clusters, 
we first need the galaxy multiplicity function. 

6.1.1. Galaxy multiplicity function 

We use the galaxy formation model described in detail in 
Valageas & Schaeffer (1999). This previous study is consis- 
tent with the present work (all parameters have the same 
values) and it was checked against observations for var- 
ious galaxy properties. In particular, it is based on the 
same formalism described in Sect. 2.1 to derive the mass 



Sx — 



L 



X 



47rr;„™(z)2 ' 



(38) 



functions of dark matter halos. We briefly recall here the 
characteristics of this model we need for our purposes. We 
define galaxies by the requirement that two constraints 
be satisfied by the underlying dark matter halo: i) a viri- 
alization condition A > Ac and ii) a cooling constraint 
tcooi < tH{zform) which statcs that the gas must have 
been able to cool within a few Hubble times at formation. 
We assume that the Hubble time at formation is given 
by the dynamical time, i.e. that the density of these dark 
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Fig. 14. The comoving galaxy B-band luminosity function 
at z = 0. The solid lines show our scaling approach and the 
dashed-lines a modified PS prescription. The dotted lines 
correspond to the scenario (C) where galaxies are defined 
by a constant density threshold Ac(z), for both prescrip- 
tions (thus the upper dotted curve is the standard PS mass 
function) . The data points are observations from Loveday 
et al. (1992) (circles), EUis et al.(1996) (filled squares) and 
Zucca et al. (1997) (triangles). The upper (resp. lower) 
panel corresponds to the SCDM (resp. OCDM) cosmol- 
ogy- 



matter halos does not evolve much after collapse. Note 
that just-collapsed halos (A = Ac) satisfy tdyn ~ tn by 
definition. Thus, the dark matter radius of the halos we 
identify with galaxies is given by the conditions: 



^cool < td yn 



and 



p > (1 + Ac)p. 



(39) 



In other words, the radius i?(T, z) of galaxies of virial tem- 
perature T at redshift z is: 



RiT,z) ^ Mm {RcoohRvr 



(40) 



where Rcooi is a "cooling radius" defined by tcooi ~ tdyn 
while Ryir is the "virial radius" defined by A = Ac. Note 
that for clusters, by definition, we have R = Rvir- At 
low redshift, for small masses, cooling is very efficient so 
that the virialization condition is the most constraining 
one in (p9). Hence R = R^r and these galactic halos are 



defined by the usual density contrast threshold Ac(2;). On 
the other hand, for large masses (i.e. high T) cooling is 
inefficient and only occursfor high gas densities, so that 
galactic halos are defined by i? = Rcooi- This means that 
their mean density contrast is larger than Ac and that 
these objects formed at a larger redshift than the one we 
would have obtained by considering just-collapsed objects 
defined by A = Ac. At large T, where bremsstrahlung 
is the main cooling process, the cooling radius goes over 
to a constant Rcooi ~ 100 kpc. Finally, we use a simple 
star formation model which takes into account infall into 
the inner parts of the galaxy, feedback from supernovae 
(proportional to 1 /T for small galaxies as in Kauffmann et 
al. 1993) and with a star formation time-scale proportional 
to the dynamical time (see Valageas & Schaeffer 1999 for 
details). 

We recall in Fig.|lJ the B-band luminosity function we 
obtain at z = 0. We can check that our predictions agree 
reasonably well with observations. Moreover, we can check 
that the "extended PS" prescription (dashed lines) pre- 
dicts more intermediate mass halos and fewer very mas- 
sive objects than the scaling mass function, as noted in 
Sect. 2.1. This is consistent with the behaviour we already 
found for clusters in Scct. ^.^ . In addition, in order to 
clearly show the importance of correctly identifying the 
galactic halos we also display in Fig.|l^ the luminosity 
functions we would obtain (dotted curves) for a "model" 
(C) where all halos are defined by the virialization con- 
straint (constant density contrast): 



(C) : Agaiix) = Ac(z) for aU x 



(41) 



with all other parameters (i.e. the star formation model) 
kept unchanged. Of course, this "model" is only shown 
for illustrative purposes since it is clearly inadequate for 
massive halos. Indeed, in the case (C) at large masses 
we identify clusters or groups and not galaxies! Thus, we 
can see in Fig.^ that for faint luminosities which corre- 
spond to small halos, the "model" (C) superposes onto 
our actual galaxy model while for large luminosities it 
leads to huge galaxies which are not observed. In partic- 
ular, it implies a bright magnitude tail of the luminosity 
function which is much too flat. Note that Monaco et al. 
(2000), u sing a model for quasars similar to ours (see be- 
low Sect. |6.2.1 ), also manage to recover observations for 
both galaxy and quasar luminosity functions. However, 
in a standard procedure that has been used (Schaeffer & 
Silk, 1988) as soon as the PS prescription became pop- 
ular, in order to correct the PS mass function so as to 
use it for galactic halos they multiply the PS prediction 
by a factor exp[— (M/Mcoo;)^^^] which is fitted to the cut- 
off of the observed galaxy luminosity function. Although 
this procedure can improve the mass functions it does not 
really deal with the "subclustering problem" itself (each 
massive halo still corresponds to a cluster). In contrast, 
our prescription has the serious advantage of taking into 
account this "subclustering problem" in a very natural 
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fashion using physical arguments based on cooling condi- 
tions, (see (p9[)), and to count the individual galactic halos 
themselves so that one can study their internal proper- 
ties. Moreover, it is interesting to note that in addition to 
a strong cutoff of the luminosity function the use of the 
proper constraints ( |39| ) leads to a non-trivial "plateau" in 
the range —21 ^ A/_b ^ — 16 which provides good agree- 
ment with observations. 



6.1.2. Galaxy X-ray emission 



O 




36 38 40 42 44 46 
log(Lx) [ erg S-' , 0.5-2 keV] 

Fig. 15. The comoving galaxy X-ray luminosity function 
in the rest-frame frequency band 0.5—2 keV for the SCDM 
case. The solid line corresponds to the scaling formula- 
tion and the dashed line to the PS prescription, for the 
"Hot" scenario. The vertical line shows the cutoff at Tf{z). 
For the "Cold" model the luminosity function vanishes at 
higher luminosities and it is equal to the "Hot" model 
prediction at fainter luminosities. 

As explained above galaxies are defined by the two con- 



straints (39). Hence at any redshift there is a character- 



istic virial temperature Tf{z) which marks the transition 
between the low temperature regime, where R = Rvir, 
and the high temperature regime, where R — Rcooi- It in- 
creases with z and at z = we have Tf ~ 10® K. In this 
approach, large objects defined by A = Ac with a virial 




36 38 40 42 44 46 
log(Lj;) [ erg s 1 . 0.5 2 keV] 

Fig. 16. The comoving galaxy X-ray luminosity function 
in the rest-frame frequency band 0.5—2 keV for the OCDM 
scenario. Same curves as in Fig.nS. 



temperature T > Tf{z) are made of several galaxies and 
correspond to groups or clusters. They can be subdivided 
into several subunits which verify the constraints (p9[). 
Of course, halos with a virial temperature which is only 
slightly higher than Tf{z), merely consist of one galaxy 
with some gas falling from its surroundings which have 
not cooled yet. However, we shall identity just-collapsed 
halos above Tf{z) as groups or clusters, which should pro- 
vide a reasonable estimate of the transition, while we shall 
call just-collapsed halos below Tf(z), galaxies. In partic- 
ular, the mass functions of groups and clusters we used 
above were set to for T < T/(z) (this also corresponds 
to the cutoff Mi{z) which appeared in (|9|) to obtain (?/)). 
However, this has no effect on the quantities we have stud- 
ied so far because we considered high virial temperatures 
T > 0.5 keV at low redshifts z < 1. Since in our model 
galaxies below and above Tf have a rather different history 
we consider each regime separately. 

First, we write the galaxy luminosity (due to the cool- 
ing of the gas) as: 



L 



E 



3 MhotkT 



bol 



icool 2 IXTTlptcool 



(42) 
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where Mhot is the mass of hot gas (i.e. of the order of 
the virial temperature T) and tcooi is its coohng time. 
Galaxies with T < Tf{z) are defined as "just-collapsed" 
objects (A = Ac) and they satisfy tcooi < tniz). Although 
the latter, by definition, have undergone a major merging 
event in their recent past (at a redshift z + 5z with 8z/{l + 
z) ^ 1), all of them are not necessarily in the midst of 
such a process. The time which has elapsed since the last 
merging, is measured by tH^ so we write for the galaxy 
X-ray luminosity function 



T <T 



f 



viLboi) 



dL 



bol 



^bol 



H 



r,{M) 



dM 
'W' 



(43) 



This ensures that we only count the galaxies where the 
gas has not had time to cool significantly since the last 
major merging event (while we neglect the contribution 
of the small objects which had time to cool). Then, the 
mass Mhot of hot gas characteristic of these objects is 
of the order of the total mass of gas Mg, hence we take 



Mho 



Mg (and the cooling time is evaluated using the 



gas density pg = {1 + Ac)pb)- On the other hand, massive 
old galaxies with T > Tf have already had time to cool 
since they satisfy tcooi = tH{zform) < tniz). However, 
if cooling is inhomogeneous and there is an approximate 
pressure equilibrium, a diffuse gaseous hot component may 
be present in the galactic halo with a density given by the 
condition tcooi — tniz), as we discussed in Sect. 5.1.3 for 
clusters. This component could be a left-over of the initial 
baryonic content of the galaxy (which gradually cools with 
time and falls into the inner parts of the galaxy to form 
stars) but it could also be replenished by the ejection of 
matter from the center of the galaxy by supernovae. In 
order to investigate the range of the galactic contribution 
to the overall X-ray emission we consider the following two 
models. First, we assume that there is no hot component 
(or it is negligible) so that in this "Cold" scenario old 
galaxies with T > Tf{z) no longer emit in X-rays: 



(Cold) : Lfoo, = for T > Tf{z). 



(44) 



This implies a sharp cutoff for the galaxy X-ray luminos- 
ity function at the luminosity which corresponds to the 
transition Tf (i.e. there is an upper bound for Li,oi)- Note 
however that these galaxies contribute to the X-ray lumi- 
nosity function at the redshift z/orm, when they formed. 
Second, we consider a "Hot" model where there is a dif- 
fuse hot component at the virial temperature T with a 
characteristic density: 



(Hot) 



Phot — 



Sp^nipkT 
2pAc{T)tH 



for T>Tf{z). 



(45) 



This gives the mass Mhot at the galactic radius R and 
the luminosity Lboi from ( ^^ where we now have tcooi = 
tniz). Then, we obtain the X-ray luminosity function from 
the galaxy multiplicity function described in Sect. |6.1.i| as: 
r]{Lboi)dLboi/Lboi — ri{M)dM/M. In practice, we can ex- 
pect the actual X-ray luminosity function in our universe 



to be somewhere in between these two ideal scenarios. In 
particular, since the cooling radius is equal to the actual 
radius R of the galaxy (by definition), there is not neces- 
sarily a large reservoir of hot gas at radii larger than Rcooi 
to ensure that the hot phase occupies all the volume of 
the galactic potential well, in contrast to the case encoun- 
tered for clusters. However, there may be some infall from 
the surrounding IGM (in addition to SNe ejecta). Hence 
the "Hot" scenario can be interpreted as an upper bound 
for the X-ray galaxy counts and the "Cold" scenario as a 
lower bound. Finally, we obtain the X-ray luminosity Lx 
in a given frequency band vi — V2 as in (^), which allows 
us to derive the galaxy luminosity function in this band. 

We show in Fig.|l^ and Fig.|l^ the galaxy X-ray lumi- 
nosity functions we obtain in this way for both cosmologies 
in the rest-frame frequency band 0.5 — 2 keV. The solid 
lines (resp. dashed lines) correspond to the scaling formu- 
lation (resp. the PS prescription) for the "Hot" scenario 
(psf). For the "Cold" model, the luminosity function van- 
ishes at high luminosities above the cutoff Tf{z) shown by 
the vertical solid line while it is equal to the "Hot" model 
prediction at fainter luminosities. As explained above the 
cutoff Tf{z) increases with z which leads to a larger lu- 
minosity cutoff Lf{z) at higher z. Below Tf the prefac- 
tor tcooi /in diminishes the contribution of small and faint 
galaxies since an increasingly large fraction of these ob- 
jects is cold. Moreover, the contribution of small galax- 
ies (which have a low virial temperature T) is strongly 
suppressed by the factors exp{—hpi'/kT) which enter the 
X-ray luminosity (^. 

As usual, the redshift evolution is smaller for the open 
universe (see Valageas & Schaeffer 1999) as the galaxy 
multiplicity function evolves more slowly (for the same 
reason as for the cluster multiplicity function). However, 
the cutoff Lxf{z) increases faster with z than for the 
SCDM case. This is due to the fact that the tempera- 
ture attached to these galaxies (in particular Tf) is some- 
what smaller for the low-density universe which implies 
that the factors exp{~hpi'/kT) are more sensitive to Tf 
and evolve faster with z. Note that the X-ray luminosity 
Lxf{z) is indeed lower for the open universe. As explained 
in Sect. 6. 1.1 the cooling constraint plays a greater role at 
small redshift, as shown by the comparison between the 
"Hot" and "Cold" models. This implies that the difference 
between both scenarios is largest at 2; = 0. Note that even 
for the "Hot" model the high luminosity cutoff is stronger 
at lower z. This is due to the gradual decline with time of 
the mass of hot gas Mhot attached to these halos together 
with the larger Hubble time t^iz), see ( p^ ) and (p2|). 

6.2. Quasar X-ray luminosity function 

In addition to galaxies, groups and clusters, quasars are 
X-ray emitters. In fact, as we shall see in the next sec- 
tion their contribution to the X-ray fiux is much larger 
than the galactic emission because of their harder radia- 
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tion spectrum (which is roughly similar to a power-law as 
opposed to the black-body spectrum of stars). Thus, it is 
of interest to estimate the quasar source counts since they 
dominate at fluxes smaller than those corresponding to 
clusters. Moreover, it allows us to obtain a complete de- 
scription of X-ray objects within the framework of a uni- 
fied model. Indeed, since in our model QSOs correspond 
to galaxies where some gas is accreted by a central black 
hole the quasar luminosity function is derived from the 
galaxy multiplicity function. Hence it is fully consistent 
with the cluster and galaxy mass functions we obtained in 
the previous sections. We now briefly describe our model 
for quasars, which is similar to the one used in Valageas 
& Silk (1999a,b) to study the reheating and reionization 
history of the universe. 

6.2.1. Quasar multiplicity function 

First, we assume that the mass Mq which is available to 
fuel the quasar is proportional to the sum of the mass of 
central gas Mgc and of stars which formed lately AMs: 



Mr 



F{Mg. 



AM, 



(46) 



Thus, at late times when most of the gas is consumed into 
stars the total mass of the black hole is Mbh ^ FMg. 
We use F = 0.005 since observations from Magorrian et 
al. (1998) find that Mbh ^ 0.005M^p?, where M^ph is the 
mass of the stellar bulge. From the model of star formation 
described in Valageas & Schaeffer (1999) we have: 



and 



1 



T, 



SN 



AM, = M, min 1 



tH dMs 
M, dt 



A e-^Mb, 



(47) 



(48) 



where Mb is the mass of baryons in the galaxy and as a 
function of our scaling parameter x 



X{x) 



P_ 



1 



TsN 



'{\ + A)gal{x) 



(1 + Ac 



(49) 



Here p/ 13d — 0.5 is a parameter which enters the definition 
of the dynamical time, while Tsn ~ 10® K describes the 
ejection of gas by supernovae and stellar winds (see also 
Kauffmann et al. 1993). The coefficient A measures the 
efficiency of star formation. Thus, small galaxies with T <C 
Tsn have A ^ 1 because supernovae eject a large fraction 
of the gas out of the inner parts of the galaxy. Hence 
Mgc < Mb and AM^ ^ Ms<.Mb as seen in (|^ and (||). 
On the other hand, massive old galaxies with a density 
contrast Agai > A^ have a small dynamical time, hence a 
small star formation time in our model. This leads to the 
factor ^(1 -|- A)gai in (^). Hence they have A 3> 1 and 
they have already consumed most of their gas, so that the 
quasar runs out of fuel, see ( |47| ) and (^8|). Next, if wc note 



/ecI the Eddington ratio we write the quasar luminosity 
Lq = /EdiEd as: 

Mo c2 



iQ = /i 



Ed- 



with = 4.4 10*yr, 



while the quasar life-time tg is: 



/Ed 



(50) 



(51) 



where eq = 0.05 is the quasar radiative efficiency. Finally, 
we write the quasar multiplicity function we would obtain 
without any scatter as: 



dMQ 

Mq 



tQ dM 



(52) 



where rjJjv'I)dM / M is the galaxy mass function obtained 



in Sect. 6.1. Here Im ^ tn is the evolution time-scale of 



galactic halos of mass M defined by: 



t~l = Max 



1 



H 



|..(M) 



(53) 



We use Aq ~ 0.05. Note that here we only have two free 
parameters (fEdF) and (Ageg/ZEd) so that we could for 
instance use a larger /Ed with a smaller F and larger 
Ag. Moreover, the properties of quasars show a signifi- 
cant scatter. For instance, Magorrian et al. (1998) find 
that the decimal logarithm of the ratio F = M bh /Msph 
obeys a Gaussian distribution of mean —2.28 and standard 
deviation ~ 0.51. Hence we assume here that the actual 
Imninosity Loscat of QSOs is related to the luminosity Lq 
defined in ( p0| ) by: 

Lg,cat-e"Lg , P{u) ^ ^L= e-^" /^''i (54) 



where P{u)du is the probability distribution of the ran- 
dom variable u. Thus, we assume P{u) to be a Gaussian 
of width cr„ = 1.15 (resp. (t„ = 1) for the SCDM (resp. 
OCDM) scenario. Note that cr„ = 1.15 corresponds to the 
scatter observed by Magorrian et al. (1998) for the ra- 
tio F. Although part of this scatter may be due to the 
observational noise, which would decrease cr„, the random 
variable u also describes the scatter of the other properties 
of quasars like the Eddington ratio /Ed- Then, we obtain 
the final quasar luminosity function from (|5l) and (^3): 



VQscat{LQs 



du P{u) TjQ {Lqscat e ") 



(55) 



Since observations suggest that the Eddington ratio 
/Ed increases with the luminosity Lq (Padovani 1989; 
Wandel 1998) we use the parameterization: 



/Ed = Min 



/ L 



Ed 



V 1047 erg/s 



0.33 



(56) 



in order to be consistent with the data obtained by 
Padovani (1989). This author finds that most of the ob- 
served dependence of /Ed on redshift is accounted for by 



Valageas & Schaeffer: Multiplicity Functions and X-ray Emission of Clusters and Groups versus Galaxies and Quasars 19 



O. 1 6<z<0.4 




H — hd — \ — — I — I — IM — M — I — I — I — — I — I — I — H 



1 .9<z<3.3 




— I — I — f — I — I — I — I — I — I — I — I — \ — — I — I — I — H 



3.5<z<4.5 





— h- 1 — \ — I — I — I ' I f } I — I — \ — \ — I — I — I — I — 



Fig. 17. The comoving quasar B-band luminosity func- 
tion for the SCDM case. The data points are observations 
from Pei (1995). The soUd hues correspond to our fiducial 
model, the dotted lines to (A), the dot-dashed Hues to (B) 
and the dashed lines to (C), see main text. 

the dependence on luminosity (low z quasars are fainter 
than high z QSOs) hence we do not add any explicit de- 
pendence on redshift in our model (p6|). However, in order 
to estimate the influence of these various parameteriza- 
tions we also consider for the B-band luminosity functions 
displayed in Fig.|l^ and Figjl^ the following two alterna- 
tive models (A) and (B) with: 



(A) : /Ed = 1 
and 

(B) : tQ = Max (2 x 10** (1 + ^)"^-^ 



(57) 



(58) 



Thus, model (A) corresponds to a constant Eddington ra- 
tio of unity while model (B) corresponds to an accretion 
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log(Lj|/Lg) 



Fig. 18. The comoving quasar B-band luminosity function 
for the OCDM case as in Fig.|l^. 



time which is roughly proportional to the Hubble time 
(with a lower bound so that /Ed < 2). Finally, in order 
to point out the importance of correctly identifying the 
galactic halos we also display in Fig.|l^ and Fig.|l^ the 
"model" (C) where all halos are defined by the virializa- 
tion constraint as in (pl[). Of course, as in Sect, 6. 1.1 we 
only show this "model" for illustrative purposes since it 
is not valid for massive halos. However, note that usual 
analytic models based on the PS approach (e.g., Haiman 
& Menou 2000) define halos as in (C). 

We show the B-band quasar luminosity functions that 
we obtain in Fig.|l^ and Fig.|l8|for both SCDM and OCDM 
cosmologies. First, we note that our result (solid lines) 
agree reasonably well with observations over the whole 
range < z < 4.5. Moreover, we can check that both mod- 
els (A) (dotted lines) and (B) (dot-dashed lines) are quite 
close to our fiducial model. Thus, although our model is 
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very crude our predictions should provide a reasonable es- 
timate of the quasar multiplicity since they are not too 
sensitive to details. Note, however, that the model (B) 
shows a slightly better agreement with observations as it 
leads to a stronger decrease of the quasar number den- 
sity at low z. This is due to the redshift evolution of the 
quasar life-time Iq which implies a smaller Eddington ra- 
tio /Ed at low z. On the other hand, we can clearly see 
that the "model" (C) fails at low redshift z < 2. Indeed, 
since it counts groups or clusters as one halo it predicts a 
significant number of huge quasars, which is inconsistent 
with observations. Of course, at high redshift where the 
cooling condition (first constraint in (|39|)) plays no role 
the "model" (C) superposes onto our model since all 
galaxies are defined by the usual density threshold Ac(z). 
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Fig. 19. The comoving quasar X-ray luminosity function 
in the frequency band 0.5 — 2 keV for the case fim = 1- The 
data points are observations from Miyaji et al. (1998). 

Thus, one cannot use the (standard) PS prescription 
to obtain the QSO luminosity function at low z. Hence one 
should be careful about the use of the PS mass function to 
derive the redshift evolution of quasar properties from the 
observed luminosity function. For instance, while Haiman 
& Menou (2000) find that this procedure implies that the 
Eddington ratio /sd or the mass ratio F must drop by a 
factor ~ 100 at low z between ^ z ^ 3 (which would 
steepen the dashed curves shown in the upper panel in 
Fig.|l^ and Fig.^, they note that a possible "solution" 
would be to correct the halo mass function in order to 
make sure that one counts galactic halos and not clusters. 
Indeed, our results show that the correct procedure (i.e. 
identifying galactic halos) provides by itself a reasonable 
agreement with observations, without any additional ad- 
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Fig. 20. The comoving quasar X-ray luminosity function 
in the frequency band 0.5 — 2 keV for the OCDM scenario. 

hoc redshift dependence. It could be tempting to introduce 
a redshift dependence of the form (|5^), to improve some- 
what the fit to observations. Nevertheless, we consider the 
present models too crude to allow one to draw such con- 
clusions from the observations: such a mouve would plug 
theoretical inaccuracies into an artificial "observed" evo- 
lution of ( |58| ) with redshift. Note however that such an 
evolution with redshift may exist. 

6.2.2. Quasar X-ray emission 

From the quasar multiplicity function we obtain the X-ray 
source counts which are identified as QSOs as in ( ^8|) . We 
assume that the quasar spectrum is locally a power-law 
Li, cx v~^^^ around v\ = \ keV, normalized by {L\l L\,oi) — 
0.028 with Li = viL^[vi). We present in Fig.^and Fig|2^ 
the comoving quasar luminosity function we obtain for 
both SCDM and OCDM cosmologies. The frequency band 
0.5—2 keV corresponds to the observed spectrum (i.e. light 
was emitted between 0.5(1 -I- z) and 2(1 -I- z) keV). We see 
that we obtain a reasonable agreement with observations 
from Miyaji et al. (1998), both at low redshift z ~ 0.3 
and high redshift z ~ 3.5. This could be expected from 
the results of Fig.|l^ and Fig.^ for the B-band luminosity 
function. 

6.3. Galaxy and quasar versus group and cluster counts 

From the X-ray luminosity functions obtained in the pre- 
vious sections we can derive the surface density on the 
sky of sources brighter than an X-ray flux limit Sx , tak- 
ing into account the contribution from galaxies as well 
as from QSOs, groups and clusters. We show our results 
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Fig. 21. Upper panel: the number of sources per square 
degree brighter than an X-ray flux hmit Sx for the case 
flm = 1- The soHd lines correspond to galaxies for the 
"Hot" model (upper curve) and the "Cold" model (lower 
curve). The dashed line corresponds to groups and clus- 
ters. The dotted line shows the quasar number counts. We 
only display the scaling model. The data points for quasars 
are from Hasinger et al. (1998) (disks for the 1112 ksec 
PSPC and triangles for the 207 ksec HRI) and McHardy 
et al. (1998) (crosses). The data points for clusters are as 
in Fig.|l^. Lower panel: same curves for an open universe 
with rim = 0.3. 



for the scaling model in Fig.pl|. First, we note that the 
prediction of the "Hot" model for the galaxy counts is 
much larger than for the "Cold" scenario, especially for 
the low-density cosmology. This could be expected from 
the comparison of Fig.|l^ with Fig.|l6|. Thus, in the "Hot" 
scenario the contribution from individual galaxies dom- 
inates over groups and clusters at Sx ^ 10~^^ erg s^^ 
cm~^ while in the "Cold" scenario this only holds for 
Sx ^ 10"^^ erg s~i cm'^ for the SCDM case and for 
even smaller fluxes for the OCDM cosmology. This is due 
to the much lower cutoff Lxf in the OCDM scenario, as 
seen in Fig.p^ and Fig.^. Of course, the characteristic 
fluxes of galaxies are much smaller than for clusters be- 
cause of their smaller mass and temperature. Moreover, 
a signiflcant part of the gas of bright galaxies has been 



able to cool due to their higher redshift of formation. This 
leads to the sharp high-flux cutoff of the galaxy counts 
as compared with the group and cluster counts, even for 
the "Hot" model. However, at low Sx galaxies provide 
a larger contribution than groups because they are much 
more numerous. 

In the "Cold" scenario, the observation of these X- 
ray emitting galaxies would provide a direct signature of 
galaxy formation and of the associated cooling process. 
Bright galaxies already start to appear at Sx ~ 10~^^ 
erg cm~^, typically one per square degree, and are 
within the AXAF sensitivity limitsQ In our model (see 
Valageas & Schaeffer 1999), such bright galaxies typically 
correspond to a virial temperature of 10^ K and a baryon 
mass of 2.3 x lO^^M©. They are just forming at z = 3, 
with an X-ray luminosity of 1.5 x 10'*'^ erg/s over a time 
of 1.4 X 10^ years. However, in the "Hot" scenario most 
X-ray emitting galaxies would have already seen a large 
part of their gas content cool and the observed luminos- 
ity would arise from a diffuse hot component, left-over 
from the last merging event or replenished by supernovae, 
which only contains a small fraction of the total baryonic 
mass. In any case, it is clear that observations (and even 
the lack of galaxy detections at Sx ^ 10"^^ erg s~^ cm^^) 
would provide very interesting informations on galaxy for- 
mation and evolution. As we noted above, however, due 
to theoretical and observational ambiguities the most rele- 
vant curve is probably the sum of both contributions from 
galaxies and groups. Thus, we see that the predictions ob- 
tained for both cosmologies are similar. 

On the other hand, because of their harder radiation 
spectrum, quasars provide an important contribution to 
the X-ray source counts, which actually dominates over 
the whole range 10~^^ < Sx < 10~^^ erg s~^ cm~^. We 
note that although we recover the right abundance of AGN 
source counts at high luminosities Sx ^ 10^^** erg s~^ 
cm~^ we somewhat underestimate the number of low lu- 
minosity objects Sx ^ 10~^^ erg cm~^. This may 
suggest that a more detailed model of QSOs is needed 
in order to match exactly the observations. On the other 
hand, we note that McHardy et al. (1998) flnd that at low 
fluxes a new population of sources appears, which consists 
of narrow emission lines galaxies which could partly corre- 
spond to starburst galaxies, which we did not speciflcally 
include in our model. 



7. Summary and Conclusion 

We have examined the predictions of the scaling model 
for the number of X-ray emitting objects as a function 
of redshift. This approach gives the multiplicity of non- 
linear structures directly from the non-linear density field 
and should be considered as an improvement over the PS 



^ We are indebted to R. Mushotzky for a discussion of this 
point. 



22 Valageas & Schaeffer: Multiplicity Functions and X-ray Emission of Clusters and Groups versus Galaxies and Quasars 



approximation. This model gives more precise counts and 
can also be used for any density contrast, both much larger 
than the usual value of 200 as well as much smaller. It also 
solves the hierarchy (cloud-in-cloud) problem and allows 
us to describe structures embedded within other virial- 
ized condensations. Thus we can simultaneously describe 
galaxies (and quasars) as well as galaxy clusters. 

For the cluster temperature and X-ray luminosity dis- 
tributions, the counts obtained in this way typically differ 
by a factor of two from the PS approximation, but the 
differences may reach an order of magnitude for extreme 
cases (very large or very small objects). With an initial 
CDM spectrum that is normalized in the same way as in 
the numerical simulations that reproduce the same data, 
the scaling model reproduces the currently available obser- 
vations, while the PS approximation docs not (the latter 
would, provided the normalization of the power spectrum 
is modified, with however a different modification depend- 
ing on the observations to be reproduced). The evolution 
with redshift is different too, for reasons that are well un- 
derstood. The counts for T > 5 keV, for instance, peak 
at z = 0.4 rather than z = 0.2 for fim = 1, with a much 
larger normalization. 

This allows a more accurate estimation of the Sunyaev- 
Zel'dovich distorsion parameter y along a given line of 
sight (summing over the contribution of clusters). For 
the same CDM initial conditions as above, we get {y) = 
2 X 10~^ that induces fluctuations {Sy"^)^/"^ ~ 10^^ in the 
CMB for = 0.3, and values which are a factor of two 
smaller in a critical density universe, the bulk of the con- 
tribution arising from clusters at z < 1. 

Then, we have shown that models of the intra-cluster 
medium which include a characteristic temperature Tad ^ 
0.4 keV reproduce the main features of the observed tem- 
perature - X-ray luminosity relation, independently of the 
details of the models. This provides a robust estimate of 
the luminosity function which is seen to show a very weak 
dependence on redshift. 

Next, we have recalled how to build a model for galax- 
ies and quasars in a way that is consistent with the de- 
scription of clusters and of the underlying density field. In 
particular, we have pointed out the advantages of our ap- 
proach which allow us to study these high-density objects 
which cannot be dealt with by the standard PS prescrip- 
tion. Thus, we show that the analysis of the observed QSO 
luminosity function in the light of the PS mass function 
leads to discrepancies and to erroneous conclusions due to 
the intrinsic limitations of this theoretical approach. On 
the other hand, thanks to its more extended range of ap- 
plications, our method can be meaningfully used to study 
these objects and it provides a much better agreement 
with the data. This also holds for the galax;y luminosity 
function, for the same physical reasons. 

Finally, we have taken advantage of being able to 
use the same basic model to draw a global description 
of the X-ray emission from all structures. Thus, we find 



that quasars dominate the X-ray source counts over the 
available range of fluxes. Clusters and groups provide a 
non-negligible contribution (^ 1/5 of the quasar counts) 
over the same range. On the other hand, at low fluxes 
Sx ~ 10"^® erg s~^ cm~^ one starts to get access to indi- 
vidual galaxies, especially for a critical density universe. 
For the OCDM scenario, this only occurs if there is inho- 
mogeneous cooling so that a hot diffuse gaseous compo- 
nent remains in the galactic halo while gradually losing 
mass in the form of cold gas (otherwise, if there is no such 
hot gas, galaxies only appear for Sx ~ 10~^* erg s~^ cm~^ 
in this OCDM cosmology). These objects, that were more 
luminous in the past and start appearing in deep enough 
surveys, are undoubtedly a new challenge. In particular, 
they should be accessible with the current sensitivity of 
AXAF. Thus, observations of X-ray emitting galaxies (or 
the lack of detection) will provide interesting information 
on galaxy evolution. They will constrain the amount of hot 
gas within galactic halos (which in turn could give some 
constraints on the feedback from supernovae and the in- 
fall from the IGM) and especially uncover some massive 
galaxies while they are just being formed. 
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